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Radiative transfer in a clumped winds is used to describe X-ray emission line profiles observed in the XMM- 
Newton RGS spectrum of the OI star £ Puppis. ft is shown that this X-ray spectrum can be explained as 
originating from a multi-temperature collisional plasma located in the wind acceleration zone. The X-rays are 
attenuated in the clumped stellar wind, which gives characteristic profiles to the emergent lines. We specifically 
study the Nvn emission line in the spectrum of CPuP- Long RGS exposures reveal that the Nvn line profile 
is structured. On the basis of our £Pup atmosphere model, we rule out the presence of Nvn in the cool wind 
component. We suggest that the detailed Nvn line structure is due to self-absorption in the hot plasma. Wind 
clumping also affects the transfer of ionizing radiation in high-mass X-ray binaries (HMXBs) . We derive analytical 
formulae for the ionisation parameter in dependence on the parameters of wind clumping. 
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Figure 1. Model X-ray emission line profiles. 
The same stellar parameters were assumed for all 
three profiles, except of different clumping prop- 
erties as indicated. 



1. Clumping of the stellar wind 

Stellar winds of O-type stars are strongly 
clumped as confirmed by the growing number of 
observational evidences. The observed variabil- 
ity in the He n 4686 A emission line in £ Pup was 
explained as an excess emission from the wind 
clumps [1]. The Ha variability in a large sam- 
ple of O stars was reproduced by a wind model 
consisting of coherent shells [2] . 

The empirical mass-loss estimates shall be gen- 




Figure 2. The EUV spectrum of £ Pup, observed 
with FUSE (thin line), compared to a PoWR 
model spectrum (dotted). The resonance dou- 
blets of Civ, Nv and Ovi are well reproduced, 
as well as the "forest" of iron-group lines. The 
O vi doublet can only be fitted with models when 
assuming that a diffuse X-ray field causes addi- 
tional ionisation. 



erally reduced by factors of a few in clumped 
winds [3]. From a study of FUV spectra of 40 
Galactic O-stars it was concluded that their mass- 
loss rates should be reduced by up to one order of 
magnitude due to wind clumping [4]. An analy- 
sis of the X-ray emission line profiles from S Ori, 
C Pup, and C Ori [5-7] has shown that the atten- 
uation by the stellar wind is significantly smaller 
than expected from homogeneous wind models. 
The observed attenuation of X-rays in the massive 
binary j 2 Vel is much weaker than expected from 
smooth stellar wind models, implying strong wind 
clumping [8]. Similar conclusions were reached 
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Figure 3. Part of XMM-NewtoriRGS spectrum of £P U P- Overplotted are two wabsxapec models with 
different temperatures. The level of continuum and the line ratios are in agreement with a plasma in 
collisional equilibrium. 



from the X-ray observations of WR 140 [9] . 

The X-ray emission line spectra of several wind- 
fed high-mass X-ray binaries (HMXB) can be ex- 
plained as originating in a wind where cool dense 
clumps are embedded in rarefied photoionizcd gas 
[10]. The stochastic variability of X-ray light- 
curve observed by XMM-Newton in 4U 1700-37 
were explained as the result of the neutron star 
feeding by a strongly clumped stellar wind [11]. 

Clumping of the stellar wind is thought to be 
a consequence of the intrinsic instability of the 
radiative-driving mechanism [12]. The ID hy- 
drodynamic modeling predicts that line-driven 
stellar winds are strongly inhomogeneous, start- 
ing from close to the photosphere [13]. Strong 
gas compression leads to the formation of cool 
clumps. The space between clumps is essentially 
void, but small gas parcels are ablated from the 
clumps and accelerated to high speed. When a 
parcel catches up and rams into the next outer 
clump, the gas is heated and emits thermal X- 
rays. The X-rays propagate trough the wind and 
suffer absorption by the cool clumped compo- 



nent of the wind. gets. Thus the wind consists 
of two disjunctive components: hot (~ 10 7 K) 
gas parcels emitting X-rays and compressed cool 
(~ 10 5 K) fragments that attenuate this radia- 
tion. 

Radiative transfer in a clumped wind is differ- 
ent from transfer in a homogeneous wind [14]. 
Formation of X-ray emission lines in clumped 
winds of massive stars was extensively studied in 
[15-17]. The main findings of these studies are 

- The opacity in a clumped wind is reduced com- 
pared to a homogeneous wind of the same mass- 
loss rate. 

- The wavelength-dependence of the opacity is 
weaker in clumped winds; if clumps are optically 
thick, the opacity becomes grey. 

- X-ray emission line profiles are nearly symmetric 
and have similar shapes across the spectrum. 

- The line profiles observed in X-ray spectra of 
single O stars can be reproduced in the framework 
of a clumped wind model. 




Figure 4. Lines observed in CP U P (overplotted are RGS1 and RGS2). The solid lines are the clumped 
wind model with M — 8.8 x 10~ 6 Af Q yr -1 . The cool-wind opacity is calculated from a PoWR model. 
The rest- frame wavelength is denoted Aq. The rest- frame wavelengths of blended lines are denoted Ai. 



Fig. 1 shows model line profiles for smooth and 
clumped winds. Two models of clumped wind 
are shown: the wind which consists of the spher- 
ical and of the flattened clumps. The line flux 
in the Fig. 1 is normalized such that it is unity 
for a line which is not affected by absorption. As 
can be seen from Fig. 1 the emergent flux in the 
line is lower in a smooth wind that in a clumped 
wind model. It is assumed that the total num- 
ber of clumps in the wind is constant. Clumps 
form at random times near some inner boundary 
and propagate radially according to a monoton- 
ically increasing velocity law. An average time 
separation between clumps is equal to the wind 
flow time Tfl = R^/v^ ~ lhr. The line profiles 
are affected by the assumed clump geometry A 
model with spherical clumps assumes that the op- 
tical depth across a clump has no angular depen- 
dence (see [18]). Slab-like clump models assume 
clumps, that are geometrically thin and aligned. 
In the latter model the optical depth across the 
clump is angular-dependent (see [15-17]) and the 
resulting line profile is more symmetric. Future 
2D hydrodynamic models will shed new light on 
the clump geometry. 



2. Stellar atmosphere model 

While our 2D clumped-wind model is used to 
calculate the formal radiative transfer, the opac- 
ity of the cool medium is obtained with the 



Potsdam Wolf-Rayet (PoWR) non-LTE code 1 . 
Shown in Fig. 2 is the observed FUSE spectrum 
of £ Pup with overplotted model. The model and 
observations are flux-calibrated. An X-ray field of 
the observed level has been included in the model 
to fit the observations. As can be seen in Fig. 2 
the P Cyg line profile of O VI is nicely reproduced 
by the model. It was the presence of such high 
ions as O vi observed in O-star spectra, which led 
Cassinelli & Olson [19] to postulate the presence 
of X-rays in stellar winds, prior to their actual 
discovery. In agreement with this classical work, 
our state-of-the-art model cannot reproduce the 
observed O VI line unless in the presence of X- 
rays. 

The scattering by the resonance Ovi line in 
the wind is needed to be taken into account for 
the analysis of the TZ ratio between forbidden and 
intcrcombination lines in helium- like Mgxi [20]. 
The transition between 2 3 5i — > 2 3 Pi levels at 
A 1034.31 A is between Ovi doublet at A 1031.91 
A and A 1037.61 A. Thus, the full radiative trans- 
fer is needed to evaluate the effect of radiation 
field on the lines ratio in helium-like ions. This is 
fully accounted for in our atmosphere models. 

From the line ratios in helium-like ions ob- 
served by Chandra and XMM-Newton in £ Pup we 
infer the X-ray formation region to be located ap- 
proximately between 1.5 R* and 10 R* [17]. This 
is consistent with the radii of formation obtained 



1 www.astro.physik.uni- 

potsdam.de / ~ wrh/PoWR/powrgridl . html 



L.M. Oskinova, W.-R. Hamann, A. Feldmeier 




28.5 



Figure 5. 
CPup. 
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Observed Nvi line in spectrum of 



independently from line profile fits. Detailed 
analysis of the forbidden-to-intercombination line 
ratios in £ Pup, £ Ori, t Ori, and S Ori is presented 
in [20]. 

3. Analysis of XMM-Newton RGS spectrum 
of C Pup 

We have analyzed the archival RGS spectra 
of C Pup. General properties of £ Pup's X-ray 
spectrum are consistent with a multi-temperature 
plasma in collisional equilibrium (CIE), as con- 
firmed by a study of differential emission measure 
distribution [21]. Fig. 3 shows part of the RGS 
spectrum with overplotted wabs*apec CIE models 
for two temperatures. As can be seen, the ratio 
of Ovm Ly/?(A16.20)/Lya(A18.99) is consistent 
with a plasma temperature of A:Tx ~ 0.6 keV. On 
the other hand, the strength of O vn helium-like 
line indicates presence of a kTx ~ 0.2 keV plasma. 
Importantly, the level of continuum seen in this 
high quality spectrum is in agreement with the 
inferred plasma temperatures. 

Continuum was only marginally detected in 
previous less sensitive observations [22] and its 
alleged weakness was used as an argument to sup- 
port the idea, that the shocks in O star winds are 
collisionless (Rauw, Pollock & Naze, these pro- 
ceedings). A recent study of collisionless damp- 
ing in stellar winds shows that in rotating mag- 
netized stars this mechanism can indeed be im- 
portant [23]. However, from both dynamic and 
energetic arguments it was demonstrated that it 
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Figure 6. Observed Nvn lines in spectrum of 
CPup (RGS1 and RGS2 are overplotted). The 
upper panel shows the sensitivity areas of RGS1 
and RGS2 across the line. The exposure time and 
the date of observation are in left and right upper 
corners. 



can operate effectively only in stars with spectral 
type B3V or later, but not in O stars [23]. 

In order to fit the observed line profiles we 
have added the observed level of continuum to 
the model lines. Fig. 4 shows the most prominent 
emission lines in the RGS spectrum. The model 
lines were obtained with our 2D clumped wind 
model [16]. The region of X-ray line formation is 
constrained from an analysis of the fir line ratios 
in He-like ions [17,20]. Therefore, the only free 
parameter in the modeling is the fragmentation 
frequency. Our model assumes time averaged sta- 
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Figure 7. Comparison between the shape of the 
optically thick Hen line observed in WR3 and 
the N vii in the X-ray spectrum of £ Pup. 



tionary mass-loss rate and the average number of 
clumps in the wind at a given moment of time is 
constant. The fragmentation frequency is an in- 
verse time interval at which moving clumps are 
crossing some arbitrary reference radius in the 
wind [17]. We prefer to use this parameter in- 
stead of porosity length [18], since the latter is 
radius-dependent whereas the former is not. 

The model lines shown in Fig. 4 are obtained 
using mass-loss estimate which does not account 
for wind clumping [24]. The fragmentation fre- 
quency for all models is 10 _4 s _1 . This corre- 
sponds to the inverse flow time, Tfl = V00/-R*, in 
the wind of £ Pup. Using more realistic mass- 
loss rates that are smaller by a factor of few 
will result in slightly higher fragmentation fre- 
quency. The Fexvn (A15.01) line is wide and 
cannot be satisfactorily fitted unless the line for- 
mation region extends further out in the wind. 
We assumed here that the line is blended with 
Fexix. The emissivity of Fe xix (A15.08) is about 
30% of Fexvn (A15.01), assuming a temperature 
kTx ~ 0.7 keV. In the most red part of the wing 
of Fexvn, the emission from Ovm (A15.17) is 
also likely to contribute. 

The Ovm (A18.97) is easily fitted with the 
present model of a clumped wind. 

Several authors have noted that the profile of 
N VII (A24.785) is different from that of other lines 
[22,25]. It was recently suggested that this line is 
blended with Nvi (A24.898) [26]. Following this 
suggestion we model Nvn line as a blend. As- 
suming ionisation equilibrium, the ratio of N VI to 
N vii is sensitive to the temperature. Fig. 5 shows 



the observed Nvi triplet line. The measured ratio 
Q = f + i/r = 1.1 ± 0.1 corresponds to a temper- 
ature T x ~ 1MK [27]. At this temperature, the 
emissivity of N vi (A24.898) is half of the emissiv- 
ity of Nvn (A24.785) (ATOMDB v.1.3.1). Using 
these values in clumped wind model, the resulting 
fit satisfactorily reproduces the N vn line width, 
as is shown in Fig. 4. 

Figure 6 shows three different XMM-Newton 
observations of the N vii line and the sensitivity 
areas of RGS1 and RGS2. The "ragged" struc- 
ture of the line profile and the dip at A24.74A is 
clearly seen in all data sets. While the depth of 
the dip is variable, its location is relatively con- 
stant. 

We employed PoWR model to study whether 
the absorption in the cool wind can cause the ob- 
served structure of the Nvn line. Nvn can be 
produced in the cool wind via Auger ionisation of 
Nv. However, our atmosphere models show that 
N VII is absent in the cool wind, even when assum- 
ing that the X-ray field is significantly stronger in 
Q Pup wind than observed. 

The line structure seen in Nvn is typical for 
the optically thick lines in the winds of hot stars. 
Fig. 7 shows the optically thick He n line in the 
spectrum of WR3, overplotted with the Nvn 
line in Q Pup's spectrum. The optical depths in 
the resonance lines of the most abundant ions 
in the hot X-ray emitting component of stellar 
winds can be larger than unity [28]. Nitrogen is 
overabundant in £ Pup, while oxygen and carbon 
are depleted. We find that the optical depth in 
Nvn line can be large enough to result in self- 
absorption. 

4. The photoionisation parameter in the 
clumped wind of HMXBs 

Clumping of the stellar wind has a profound 
effect on X-ray emission from HMXBs [10]. We 
studied the influence of clumping on photoionised 
regions. In the following we consider the case of 
optically thick clumps. The same notations as in 
[29] are used. We assume that the wind consists of 
two components: dense clumps and a rarefied in- 
terclump medium. Applying the continuity equa- 
tion we define the mean particle number density 
as 

™( r *) = ~a ~~7 v3 ' W 

47r^mHf(r* )r£ 
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where \x is mean molecular weight and tor is the 
atomic mass unit. If the density enhancement in 
clumps is D = <ra 2 >/<n> 2 , and the volume fill- 
ing factor of clumps is fy then the particle num- 
ber density of the interclump gas is 



n ic = (1 - Df v )n 



(2) 



Assume for simplicity that the interclump gas is 
optically thin and the velocity is constant. Cor- 
recting the mean radiation intensity for the ab- 
sorption in the clumps 



Jv 



47rr 2 



f(kT) 



-r(r.) 



(3) 



where f(kT) is a function of atomic parameters 
and temperature, and f is an effective optical 
depth. When all clumps are opaque (r c i ump 3> 1), 
the effective optical depth reflects geometrical dis- 
tribution of the clumps, and is found to be grey 
[15,16] 



T = 



^07"* 



(if T c l ump > 1), 



(4) 



where hq is the fragmentation frequency. Using 
the ionisation balance equation, the ionisation pa- 
rameter in the interclump medium is 



(n)r* 1-Dfv 



(5) 



For a very high fragmentation frequency, which 
recovers the case of a smooth wind Dfy « 0, Eq. 5 
reduces to the usual definition of £. Parameter q 
[29] in a clumped wind can be written as 



K - 



(6) 



We conclude that clumping affects both the 
ionisation balance and the size of the photoion- 
ized region. A higher ionisation and a larger size 
of the photoionized region can be expected. 
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